Critical behavior of thermopower and conductivity at the metal-insulator transition in 

high-mobility Si-MOSFET's 
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This letter reports thermopower and conductivity measurements through the metal-insulator 
transition for 2-dimensional electron gases in high mobility Si-MOSFET's. At low temperatures 
both thermopower and conductivity show critical behavior as a function of electron density which is 
very similar to that expected for an Anderson transition. In particular, when approaching the critical 
density from the metallic side the diffusion thermopower appears to diverge and the conductivity 
vanishes. On the insulating side the thermopower shows an upturn with decreasing temperature. 

PACS: 71.30.+h, 73.40.-c 



Scaling theory of non-interacting, disordered, electron 
gases predicts that no rinetal-insulator transition (MIT) 
occurs in 2 dimensionsErH as temperature T — > 0. Nev- 
ertheless, what appears to be a MIT has been observed 
(at finite, though low T), first in n-Si-MOSFET's^ and 
more recently in many other 2-dimensional (2D) hole and 
electron gasesO. In the particular case of Si-MOSFET's, 
the transition is most clearly visible in high-mobility sam- 
ples, roughly p > lm 2 /Vs. As the density, n, is varied, 
there is a particular value, no, above or below which the 
resistivity p shows metallic or insulating temperature de- 
pendence respectively. For the present purposes we will 
use as a working definition that negative dp/ dT indicates 
an 'insulator', and positive dp/dT at the lowest temper- 
atures we can reach corresponds to a 'metal' (possible 
deviations from this definition and the consequences will 
be mentioned later). At n not too close to no, metal- 
lic behaviour is visible over a wide range of T, roughly 
T < Q.bEp/kB where Ep is the Fermi energy. The de- 
crease of p in the metallic state for high mobility samples 
is typically two orders of magnitude larger than can be 
accounted for by electron-phonon scattering. 

Most previous work on these systems has focused on p, 
though measurements of the compressibilityQ have also 
appeared recently. The present paper presents exper- 
imental data on the low temperature thermopower, S, 
and conductivity, a = l/p, both of which are found to 
exhibit critical behavior around tiq. Earlier, a scaling be- 
haviour was describecD for the temperature dependence of 
p{T) over a temperature range ~ (0.05 — 0.3)Ep /ks- In 
contrast, we report a different type of critical behavior for 
a. When we extrapolate our data on er, typically taken 
over the range 0.3 — 4.2 K, to the T — > limit, we find a 
power-law critical behaviour as a function of (n/riQ — 1) on 
the metallic side. In addition, at our lowest temperature 
of around 0.3 K where diffusion thermopower dominates, 
S/T appears to diverge when approaching hq. At no there 
is an abrupt change in behavior of S/T, with lower den- 



sities showing an upturn in S as T is decreased. Similar 
characterisi|ics,have long been predicted for an Anderson 
MIT in 3DCniy but such a transition should not occur in 
2D. 

The main sample used for the present p and S mea- 
surements (Sample, 1) is the same as that described in 
a previous papered and the techniques used to mea- 
sure S can also be found there. This sample has no = 
1.01 x 10 15 m -2 (as defined as above) and a peak mo- 
bility p = 1.75m 2 /Vs at T = 1.1 K. S and p have 
been measured as a function of T, down to about 0.3 K, 
at many different values of n. We have also analyzed 
independent p(T, n) data for two other samples over 
the same range of T, Sample 2 from the— same wafer 
with no = 0.99 x 10 15 m~ 2 , and Sample 313 with peak 
p = 3.6m 2 /Vs and n = 0.956 x 10 15 m -2 . The major 
experimental problem was that of measuring thermoelec- 
tric voltages with the sample in the insulating state. For 
this purpose an amplifier with input bias current < 1 pA 
and input impedance > 10 12 ri was used. With some aver- 
aging it had a resolution of 0.1 pV for source impedances 
of less than a few hundred kf2 rising to about 1 pV at 10- 
20Mf2, roughly the highest sample impedance in these 
measurements. With the sample in the metallic state, 
a Keithley 182 digital voltmeter usually gave the best 
compromise of input bias current, input impedance and 
noise. All connections to the sample had isolation re- 
sistance > 50 Gf2 and all leads were well shielded and 
filtered against rf interference. 

In the metallic region n is a linear function of gate volt- 
age and it is believed to follow approximately the same 
dependence in the insulating regionta, at least close to 
no. The results on the temperature dependence of p on 
both the insulating and metallic side are not steam, but 
are very similar to those seen in previous workoll3iU. In 
the metallic regime we have fitted our data on typically 



over the rang e 0.3 K to 4.2 K, to the equatio: 
P = Po + Pi exp (-(T /Ty ) 
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where po, pi, Tq and p are fitting constants, in order to 
evah 



B 




0.1 ' « L 1 1 1 — ■ 1 

0.8 1 2 3 4 5 6 8 

n (10 15 m" 2 ) 



FIG. 1. Density dependence of the conductivity 
in the T — > limit, of Samples 2 (solid symbols) and 
3 (open symbols). For sample 3 the data have been 
shifted by adding 1 x 10 15 m~ 2 to n. 

Figure 1 shows the results on Co = 1/po as a function 
of n. All samples follow the critical behaviour 

(n V 
a-Q = cr TO +(T S [ 1 . (2) 

The solid lines are the best fits with the following pa- 
rameters (with er in units of e 2 /h). For sample 2, 
a m = 0.2 ± 0.3, a s = 13.6 ± 0.7 and v = 0.83 ± 0.03. 
Sample 1 has an identical behaviour within experimen- 
tal error. The higher-mobility Sample 3 also follows the 
same equation with er m = 0.36 ± 0.15, a s = 34 ± 5 and 
v = 1.39 ± 0.05. These results suggest v increases with 
peak mobility but clearly more data on a variety of sam- 
ples are required. The values of a m for Sample 1 and 2 
are consistent with zero within experimental uncertainty. 
For Sample 3, a m may be finite. However, if no is al- 
lowed to decrease from 0.956 to about 0.925 x 10 15 m~ 2 , 
a fit which is indistinguishable over the range of the data 
can also be obtained with a rn = ± 0.15, <r s = 32 ± 5 
and v — 1.48 ± 0.05. A small discrepancy in no could 
arise from the identification of the critical density for the 
MIT with that density, uq, where dp/ dT changes-sign, a 
procedure which has no firm physical foundation!!!! The 
critical behavior described by Eq. (Q) with a m = is 
formally the same as that expected for a (continuous) 
Anderson transition with a mobility edge at no, whereas 
a finite a m would correspond to a (discontinuous) Mott- 
Anderson transition; neither transition should arise in 
a non-interacting 2D gasao. The inclusion of interac- 
tions along with disorder is a much more complex and 



ongoing theoretical problem (e.g. see Refs. |,| and ref- 
erences therein) and it is not yet clear if such transitions 
become possible under these conditions. Similar criti- 
cal behaviour, usually with a m consistent with zerOphas 
been seen in many 3D systems, typically with valuest2l of 
v in the range 0.5 — 1.3. There are only two previously 
reported cases related to 2D. Hanein et alx3 have made 
a similar analysis to the one above for a 2D hole gas in 
GaAs and found a linear relation between cr and n, but 
with a finite a m . Feng et al£3 have also found scaling be- 
haviour in a Si-MOSFET but it appears to be unrelated 
to J - l — L u — 
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FIG. 2. The thermopower S for Sample 1 at vari- 
ous fixed electron densities n (in units of 10 15 m" 2 ). 

We now turn to the thermopower data. A selection of 
data on S is shown in Fig. 2. At n — 8.5 x 10 15 m~ 2 
diffusion thermopower, S , is almost zero and one sees 
only phonon drag, S 9 , which pSiaries approximately as 
T 6 at the lowest temperatures.Eil As n decreases, S be- 
gins to show two distinct regions with different T depen- 
dences. At T > IK there is a relatively rapid increase 
of S, roughly cx T 3 which-4s that expected for S 9 at in- 
termediate temperatures.liil At T < 1 K, S has a much 
weaker, approximately linear T dependence indicative of 
S d becoming dominant; for n < no, this low-T behaviour, 
which is characteristic of ordinary metals, is replaced by 
an upturn in S. Concentrating on the metallic region, 
the data at lowest T are taken to give the best estimate 
of diffusion S d = aT. Fig. 3 shows that a as a function 
of n appears to diverge as n — > no- One would expect a 
divergence when Ep approaches a gap in the DOS but 
the present results-are inconsistent with this explanation 
because Hall datali3 show that in the vicinity of no the 
mobile carrier density equals n within 10%. However, 
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Eq. (0) also implies a divergence of S d . Thus, with the 
assumption of a constant density of states (DOS) Eq. (||) 
is consistent with 



a(E F ) 



E r 



(3) 



Again, with a m = this is formally equivalent 
to an Anderson transition with E c being the mo- 
bility edge. The use of the Mott relation S d = 
-{■K 2 k%T /^e){d\D^E) EF with Eq. (|) and taking 
<j m = then gives 



S" 



3e(E F - E c ) 



(4) 



This result is valid only if (E F - EAlh B T > 1; in the 
opposite limit S d tends to a,constant£ilo (~ 228 /iV/K in 
3D). Numerical calculations^ show that the approxima- 
tion of Eq. (||) gives a magnitude roughly a factor of 2 too 
large when (Ep — E c )/kBT « 2 which, for our samples, 
corresponds to A = (n — nn)/nn « 0.11 at T — 0.4 K (us- 
ing the ideal DOS. an. with an effective mass of 0.19 mnh 
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FIG. 3. Density dependence of a si S d /T 
for Sample 1. The single open symbol is for 
n = 0.97 x 10 15 m~ 2 which is just below n but k B T 
broadening should make this indistinguishable from 
no. The line is the best fit to Eq. (J^) for data with 



To simulate this saturation we add A in the denomi- 
nator (but allow it to be a variable when determining the 
best fit to the data) and, rewriting Eq. (||) in terms of n, 
we have 



S d /T 



- C / x A*+(--ir 

n 



where C — vir 2 k B / \3eE c ) is a constant expected to be 
about 32/iV/K 2 for Sample 1, again using g . If a m is 
finite in Eq. (||) , then the Mott relation shows that it will 
contribute to the denominator of Eq. also softening 
the divergence at n = uq. However, the experimental a m 
is so small that this is negligible compared to the finite-T 
effect considered here. The best fit of the data to Eq. (|5|) 
gives A = 0.15 ± 0.01 and C = (9.5 ± 1.5) /iV/K 2 and is 
shown as the solid line in Fig. 3. (As with a, the fit can 
be improved if n is slightly decreased). A is consistent 
with that expected from the argument above, but C is 
too small by a factor of about 3. However, we emphasize 
that we are comparing our results for a 2D system with 
a theoretical model of an Anderson MIT valid for non- 
interacting electrons in 3D. Some progress has been made 
on calculating.^ with the inclusion of weak interactions 
and disordered. Corrections are found which are loga- 
rithmic in T and difficult to detect in thermopower; we 
are unable to explain the observed strong density depen- 
dence in terms of the calculations. We should mention 
that we can also represent the data over the same range 
using the simple expression a — — 56/n 25 ^V/K 2 , with 
n in units of 10 15 m~ 2 , but this has no obvious physi- 
cal explanation; in particular, it does not have the form 
that we would expect for S d approaching a band edge at 
n = 0, i.e. S d oc 1/n. 

The data in the insulating regime also show a critical 
behaviour qualitatively consistent with a mobility edge. 
Thus the observed upturn of S d is expected for activated 
conduction across a mobility gap with (E c — Ep) > k B T~ 
Under these conditions the 3D Anderson model predicts!!! 
(see also the numerical calculations in Ref.Ej) 



S" 



ks 
e 



A 



E c — Ep 
k B T 



(6) 



(5) 



where A is a constant of order unity. For n < 0.79 x 
10 15 m -2 the observed minima in S d occur at T m consis- 
tent with {E c — Ep)/kBT m ~ 2. Eq. (|^) would then imply 
that the values of S d at these points should all have about 
the same magnitude. However, the observed S will have 
other contributions. In particular there will be S 9 (see 
below) and also a contribution to S d from variable range 
hopping (VRH) through localized states. (When two or 
more conduction mechanisms are present, the appropri- 
ate S d are weighted by their contributions to a). The 
T dependence of our p data and other previously pub- 
lished dataE-l in the insulating region are consistent with 
Efros-Shklovskii VRH across a soft Coulomb gap. For 
this-.pjechanism one expects S d to be a constant given 
byBE3 S d = -(k B /e)(k B T /C){dlng(E)/dE) EF where 
To can be obtained from the temperature dependence of 
p, e.g. Ref. |lj, g(E) is the background DOS, and C 
a constant w 6. If we take (d\ag(E) / dE) Ep ~ l/E F 
(implying that Ep may be in the tail of the DOS) and 
again using go to estimate Ep, we find that the calcu- 
lated S d are typically a factor of two smaller than the 
values of S observed at T m . The argumentis-pot sig- 
nificantly changed if Mott VRH is assumed .liaEJ In this 
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case 

S d a T l/3 but the 

magnitudes calculated for S are 

similar. 

As far as we are aware, the only previous work which 
attempted to follow S d into the region oL2D electron lo- 
calization was that of Burns and ChaikintH on thin films 
of Pd and PdAu. They found an upturn of S d in the 
strong localization region but no divergence at higher 
conductivities. The authors attributed their results to 
the opening o£_a Mott-Hubbard gap. In 3D, Lauinger 
and BaumannEj observed critical behaviour of a and a 
divergence of S d for metallic AuSb films, but the magni- 
tude of the latter was 2 orders of magnitude-smaller than 
seen here. Other 3D experiments on SiPcJ and NbSiEJ 
saw no divergence on the metallic side. 

For completeness, we make a few comments about 5* at 
higher T where S 9 is dominant. Little is known about the 
behavior of S 9 near a MIT but it should be present-op the 
metallic side though its precise form is not knownEj. On 
the other hand, S 9 requires conservation of crystal mo- 
mentum for electron-phopiaiL. scattering so that S 9 = 
for conduction via VRH.Effij Thus, S 9 should only ex- 
ist on the insulator side when excitation to delocalized 
states occurs. Our data show that at any fixed T > 2 K, 
S rises as n decreases but crosses no smoothly, i.e., we 
no longer see divergent behavior of S at uq. These facts 
show that activated conduction must be present for all 
densities n < no that we have investigatecL-even though 
p data (both our own and those of othersEJ) appear to 
follow the Efros-Shklovskii VRH model. 

In summary, the behavior of a and S d in the 'metal- 
lic' and 'insulating' phases in Si-MOSFETs are surpris- 
ingly consistent with a 3D Anderson MIT, though such 
a transition is not expected to occur in a 2D electron 
gas. Nevertheless, it is important to remember that to 
reliably identify the observed critical behaviour with a 
MIT requires data in the zero T limit. Although our 
analysis is based on an extrapolation to zero T, the ac- 
tual data extend only to 0.3K. Thus we should be care- 
ful not to conclude that a mobility edge or a MIT has 
necessarily been observed. Even so, the present results 
provide new information on these systems that further 
constrains any theoretical model proposed to explain the 
MIT, whether such a transition be an apparent or real 
property as T — > 0. 
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